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Abstract: We study the transition probabilities for the totally asymmetric simple 
exclusion process (TASEP) on the infinite integer lattice with a finite, but arbitrary 
number of first and second class particles. Using the Bethe ansatz we present an 
explicit expression of these quantities in terms of the Bethe wave function. In a 
next step it is proved rigorously that this expression can be written in a compact 
determinantal form for the case where the order of the first and second class particles 
does not change in time. An independent geometrical approach provides insight into 
these results and enables us to generalize the determinantal solution to the multi- 
class TASEP. 

I. INTRODUCTION 

In this paper we derive finite-time transition probabilities for the one- dimensional totally 
asymmetric simple exclusion process (TASEP) with several classes (species) of particles 
In the usual single-species TASEP in one dimension particles jump on the integer 
lattice Z independently after an exponentially distributed random time with parameter 1 to 
their right nearest neighbor site, provided that the target site is empty. Otherwise the jump 
attempt is rejected. This exclusion principle guarantees that each site is always occupied 
by at most one particle. An instantaneous configuration of this system with particles 
can therefore be represented by an ordered set of integer coordinates {x} = {xi, . . . ,xn} 
where Xj+i > Xj. An intuitive form to represent this lattice gas dynamics consists in writing 
allowed local transitions as 

AO OA (1) 

where A represents a particle and represents a vacant site. Mathematically, the dynamics 
of the TASEP can be defined through a master equation for the probability P(x, t) to find 
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a configuration x at time t. The master equation reads 

^P(x, t)=Y, K,x'P(x', t) - ti;x',xP(x, t)] , (2) 

where Wx',x is the transition rate (0 or 1 for the TASEP) to go from a configuration x to a 
configuration x'. By integrating ([2]) one obtains the solution of the master equation for any 
given initial configuration y, i.e., the conditional probability to find a particle configuration 
X at time t, given that the process started from configuration y. 

This process is a discretised version of the noisy Burgers equation js, ^ or, equivalently. 



the noisy Kardar-Parisi- Zhang equation p, |6|]. On coarse-grained scale the Burgers equation 
exhibits shock solutions and it is of great interest to understand the microscopic structure 
of the shocks on the scale of the lattice. To this end, one introduces so-called second-class 
particles jsl. We denote the first class particles by the symbol A and the second class 
particles by the symbol B. According to the exclusion principle a lattice site can either be 
occupied by an A particle, or a -B particle, or it can remain vacant (i.e. occupied by a hole, 
denoted as 0). The exchange rules are 

AO -> OA 

fiO ^ 05 (3) 
AB BA. 

All processes occur with rate unity. Observe that from the perspective of the first-class 
particles the second-class particles are not distinguishable from holes. The stationary distri- 
bution of this two-component TASEP can be constructed using pair of independent product 
measures |7i]. It is well known that a single second-class particle in a system of first-class 
particles travels with the characteristic speed of the Burgers equation. Hence it can be used 
to define a microscopic position of the shock and to study its microscopic structure js-lS]. 



It was also shown that a finite number of second class partic 



es, in a uniform stationary 
background of first class particles, form a weak bound state . The TASEP with second- 
class particles and more general two-component particle systems are also interesting from 
the point of view of non-equilibrium phase transitions [15|4l7|. 

More generally, one may consider a hierarchical multi-class particle system where for 
a particle of class k all particles of higher class behave like holes. This has interesting 
applications e.g. in the microscopic study of the so-called step initial condition where in an 
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infinite system all sites x < are occupied while the positive half lattice x > is vacant. 
This initial condition has been studied in great detail for the single-species TASEP, see e.g. 
the work by Johansson 18| for a seminal study of current fluctuations. Recently, interesting 
passing probabilities for the multi-class exclusion process have been obtained 19|. The 
stationary distribution of the multi-class TASEP has been constructed in 20| by extending 
the approach of [7] . A matrix product formalism for the stationary measure of a multi-class 



2l|. 



TASEP was discussed in 

A determinantal approach to solve 
infinite lattice has been introduced in 



or the dynamics of the single-species TASEP on the 



22| . It is based on the explicit solution of the master 



equation for the conditional probability of finding the system in a certain configuration at 
time t for a given initial configuration, using the technique of the Bethe Ansatz, originally 
employed for studying the relaxation spectrum of a periodic chain j23|. Moreover, it was 
shown that the solution of master equation for N particles on an infinite lattice can be 
written as the determinant of an x matrix. This determinantal form has allowed for 
a rederivation and extension of the work of Johansson on the step initial condition 



26| and has inspired a substantial body of further work on other dynamical properties of 



the exclusion process, see e.g. 



27H34|. In the present work we obtain an exact solution 



of the master equation for a finite number of first class and second class particles on an 



infinite lattice with arbitrary initial condition. We use the nested Bethe ansatz |35l-l37l|. thus 
generalizing the approach of 22| for the single-species case. In a recent study the technique of 
nested Bethe ansatz was used to investigate the spectral structure of the transition matrix 
for a multi-species exclusion process [38!|. In particular, the relaxation time-scale of the 
system was studied and it was shown that the dynamical exponent remains same as in the 
single-species case. In this paper, we use nested Bethe ansatz to solve the master equation 
of the two-species case. We obtain a compact determinantal representation of the solution 
for the case where first and second-class particles do not change their ordering. We also 
discuss the extension of our results to multi-species exclusion processes. 

The paper is organized as follows. In Sec. |Tl]we show how the master equation can be 
solved exactly for the two-species case using the nested Bethe ansatz. In section IIIII we 
give a formal proof of the determinantal representation for transition probabilities without 
interchange of particles. In section IIVI we present an alternative diagrammatic approach 
and discuss how using this approach the determinantal representation can be understood. 
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We conclude this paper by discussing the possible generalization to multi-species systems in 
section |Vl 



II. EXACT SOLUTION OF THE MASTER EQUATION USING BETHE 

ANSATZ 

In this section we present an exact solution of the master equation for the TASEP with 
first and second class particles. To illustrate the approach we first consider the case of total 
particle number N = 2. Each of these particles can be either of type A or of type B. 
Therefore there are 4 different particle sequences: {AA, AB, BA, BB}. 

Let Xi,X2 with 2:2 > Xi + 1 be the positions of the two particles at time t and yi,y2 
be their initial positions with y2 > Hi + 1- The quantity P^'^'^ {xi, X2;t\yi,y2;0) denotes 
the probability to find the system in sequence Q' G {AA, AB, BA, BB} with occupancies at 
xi,X2 at time t, given that at t = the sequence was Q G {AA, AB, BA, BB} with occupied 
sites at yi,y2- Note that if the initial sequence is AA,BA or BB, then no interchange 
between A and B particles is possible as the system evolves in time. For these cases the 
sequence remains the same for all times and the particles execute usual TASEP dynamics 



with transition probabilities obtained in [22|]. However, if the initial sequence is AB, then 
as a result of interchange between the A and B particles the sequence may change and the 
possible sequences at time t are AB, BA. 

In the absence of nearest neighbor pairs, i.e. for X2 > Xi + 1, the master equation becomes 

^pQ'\Q(xi,X2; t\y,, 1/2; 0) = P«'l«(a;i - 1, X2; t\y^,y2] 0) + P«'l«(xi, X2-I; t\y^, y2\ 0) 



dt 



-2P«'IQ(xi,a;2;tbi,i/2;0). (4) 



For X2 = Xi + 1, the form of the equation is different depending on Q and Q'. For Q G 
{AA, BA, BB} we have Q' = Q and the equation takes the form 

^P«l«(xi, t\y,, t/2; 0) = P«l«(xi-1, t\y,, ^2; 0)-P«l«(xi, Xi+l; t\y^, y2; 0). (5) 

dt 

For Q = AB and Q' = AB we have 

^P^^l^^(xi, t\yu 7/2; 0) = P^^l^^(xi-1, t\yu y2, 0)-2P^^I^^(a;i, t\yi, y2, 0) 

(6) 
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while for Q = AB and Q' = BA 



^P^^l^^(xi, xi + 1; ys; 0) = P^^l^^(xi - 1, a;i + 1; 1/2; 0) + P^^l^^(xi, xi + 1; t\y^, y^, 0) 



(8) 



-P^^l^''(xi,xi + l;t|i/i,y2;0). 

To solve the above set of equations for a given initial sequence Q and initial positions 
{yi, 1/2} we define 

j p^^lQ(a;,,x2;t|i/i,2/2;0)\ 
P^^l«(xi,x2;t|yi,y2;0) 
P^^l«(xi,x2;t|i/i,y2;0) 
VP^^I<3(a;i,X2;t|yi,2/2;0)y 
Notice that in order to keep the notation simple we do not mark the dependence of this 
vector on the initial sequence Q. Now we use Bethe Ansatz to write \P{t)) in the following 
form 

(27r) Jo Jo 

where 

V Sbb,q I 

The elements of the 4x4 matrix S(pi,]92) and the quantity e (^1,^2) are yet undetermined 
functions of pi and p2- Substituting this ansatz into (j4]) we find 



e(Pi,P2) 



2 - e"*P^ - e-'P^ 



(11) 



and to satisfy the boundary conditions for the nearest neighbor pairs resulting from Eqs 
([n]), dZ]), one needs 



S(P1,P2) 
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(12) 



The complete solution can then be written in terms of the functions 



gi(p+iO)j"- 



(13) 



(7) 
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introduced in 22|. From now on we will drop t from the argument of the above function 
and write it simply as Fn{x), for notational simplicity. 

For the present case of two particles the values of n that occur are 0, ±1. For Q = 
AA,BA, AB we have 

P«l« (xi, xa; t\yi, y^, 0) = Fo{x, - yi)Fo{x2 - 2/2) - Fi(x2 - yi)F_,{xi - y^) (14) 

and all other elements of \P{t)) with Q' ^ Q are zero. For Q = AB we have 

/ \ / \ 

^o(a;i - yi)Foix2 - t/2) - ^0(3^2 - yi)Fo{xi - 7/2) 
Fo{x2 - yi)Fo{xi - 1/2) - Fi{x2 - yi)F_i{xi - 1/2) 

V y V J 

(15) 

The important point to note for > 3 is that there are no new constraints from the 
boundary condition when more than two particles are on adjacent sites. This is expected 
from the integrability of the model 39|, |40| , but we have also verified this explicitly in our 
calculation (see appendix). For = 3 there are 8 possible sequences of A and B, ordered 
in a vector \P{t)) according to the following formal tensor product: 

/ aaa\ 

AAB 



P^^\^B (xi,X2;t||/i,y2,0) 
P^^l^^ (xi,X2;t|2/i, 1/2,0) 





ABA 
ABB 
BAA 
BAB 
BBA 
BBB 



(16) 



The Bethe Ansatz solution has the form 

\P{t)) ^ 



(27r)3 



2tt 



dpidp2dp3e~ 



+S?2(pi,P2)e*(^^^^+^^'^^+^^'^^^ + T.ls{p2,P3)e'^''''''^"'^'^^^'''^ 
+T.Up,,p,)T.Up,,p2)e'^^'^^^'^^^^^^ 

+E^3(pi,P2)S?2(pi,P3)SL(P2,P3)e*(^^^^"'^^^^+^^^^^}|P(0)) 



(17) 
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where T^l^iPj-iPk) = ^{Pj,Pk) ® ^2 and T.l^{pj,pk) =12® ^{Pj-,Pk)- Here I2 is the 2 x 2 
unit matrix. Here the upper index 3 refers to the number of particles = 3. One crucial 
requirement for the integrability of the model is that these S-matrices must satisfy Yang- 
Baxter criterion (see appendix for details). 

Using the definition in (JT3l) the above solution can be written in terms of Fn{xi — yj) 
with n = 0,±1,±2. For example, for the initial sequence Q = AAB we find (with x = 
(xi,X2,a;3), y = (yi, 2/2, 2/3)) 



P 



AABIAAB 



x;t|y;0) 



^o(a;i - yi) F_i{xi - 1/2) F_i{xi - ys) 
Fi{x2 - yi) Fq{x2 - 2/2) Fq{x2 - ys) 
Fi{x3-yi) Fo{x3-y2) Fq{x3 - y^) 



(18) 



P 



ABAIAAB 



x; t\y; 0) = Fo{xi - yi) {Fq{x2 - y3)Fo{x3 - 7/2) - Fi{x3 - y2)F-i{x2 - 2/3)} 
+F_i(xi - 1/2) {^^2(3:3 - ?/i)F_i(x2 - 2/3) - Fq{x2 - y3)Fi{x3 - yi)} 
+Fq{x2 - yi) {Fi{x3 - y2)F_i{xi - ^3) - Fo{xi - ?/3)^o(a;3 - 2/2)} 
+F_i(x2 - 2/2) {Fi{x3 - yi)Fo{xi - ^3) - i"2(a;3 - 2/i)F-i(xi - 2/3)119) 



pBAAlAAB ^ ^^^^^ _ y^-^ 1^^^^^ _ _ ^3) _ Fo{x2 - 2/l)F_i(xi - I/3)} 

+Fq{x3 - 2/2) {i^o(a;i - y3)Foix2 - yi) - Fi{x2 - 2/i)F_i(xi - ?/3)} 
+F2(x3 - yi) {Fq{x2 - 2/2)^^-2(^1 - 2/3) - i^-i(a;2 - 2/2)^^-1(3^1 - 2/3)} 
+Fi(x3 - 2/1) {i^-i(a;2 - 2/2)i^o(a;i - 2/3) - ^^0(2:2 - 2/2)i^-i(a;i - 2/3)^20) 

and all other components of \P{t)) are zero. Thus, starting from any given initial sequence 
Q, all components of \P{t)) can be determined. 

Similarly, for = 4, the possible 16 sequences can be represented in the tensorial ordering 




. The Bethe Ansatz solution can be written down in terms of the matrices 



^{2 {Pu Pj) = S ipupj) ® I2 ® I2 (21) 
^23 (Pi .P3) = I2 ® S {pi ,pj)^ h (22) 
(Pi, Pj) = I2 8) I2 8) S {pi,pj) (23) 

and in this way one can construct \P{t)) for any value of A^. The total of number of sequences 
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is 2^ and these sequences are ordered as in 

A ■ 

B 

The 2^ dimensional E-matrices are constructed from the tensor product of Ti{pi,pj) and I2 



matrices: 



etc. and in general 



^i2iPuPj) = ^(PhPj) ® I2 ® I2 ® ® I2 (24) 

^23{Pi,Pj) = I2 ® ^{Pi,Pj) ®h® ® I2 (25) 



wiP^^P,) = ihf'-' ® m,P,) ® (12)^"'^^-' (26) 



The conditional probability P^'^^ (xi, X2, xat; t\yi, 1/2, I/tv; 0) has the Bethe ansatz form 

I^W) = tAv r'^c^Pic^P2...rfp^ Vp(gig2...giv)e<p«i^^+^'?^^^+-+^'5^^-) \pm . (27) 

In this Bethe representation {gig2---giv} is one particular permutation of {12. ..A^} and 
'P{QiQ2---Qn) is constructed from the product of {Pk,Pj) matrices given by Eq. fl26l) . 
This product corresponds to the series of elementary permutations performed on the mo- 
menta {pi,p2, ...,Pn} to arrange them in the order {pqi,Pq2, ■■■,PQn}- The summation is 
over all possible A^! permutations. One can show by adapting the techniques of 40|] that 
Eq. ( 1271) solves the A^-particle master equation for the TASEP with second class particles. 

Notice that the matrix E in Eq. (JT2l) may be decomposed as E = D + A where D is the 
diagonal part and A the off-diagonal part which arises from the interchange AB — t- BA of 
the two species of particles. Correspondingly the product V{QiQ2---Qn) can be expanded in 
terms containing an increasing number of A matrices which represent the number of inter- 
changes in the final configuration. The purely diagonal part in this expansion corresponds to 
configurations in which the original order of first and second-class particles has not changed. 
We investigate these realizations of the stochastic dynamics in more detail in the following 
section. 



III. THEOREM FOR NO INTERCHANGE 



In the previous section we have shown that using Bethe Ansatz it is possible to solve the 
master equation for any number of particles with any sequence. The Bethe Ansatz form 
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of the solution for particles has, in general, A^! number of terms. In the single-species 



TASEP case, it was shown in 



22[ | that the solution of the master equation can be written 



as a determinant of a matrix F whose elements are Fij = Fi-j{xi — yj) with Fi^j{xi — yj) 
defined according to ( fT3l) . In this section we show that there is compact determinant 
representation of the transition probability for the case where no interchange between A 
and B particles takes place, i.e., where the sequence at time t is the same as the initial 
sequence at time t = 0. We state the main result: 



Theorem for no interchange: Consider a system of A^ particles, each of species A or B, 
starting with an initial sequence Q at time t = taken from the tensorial ordering 




The probability P<3IQ 

(xi, X2, xat; ?/2, 1/Af; 0) where the initial and final sequence 
remains the same can be written as the determinant of the N x N matrix G, where 

Gij = Fn(i,j) {xi - yj) (28) 

with matrix elements F„ defined in (fT3ll and n{i,j) = sgn{i — j) {\i — j\ — uab) with uab 
being the number of AB pairs occurring between the z-th and j-th particle of the sequence Q. 



Remark: To illustrate this result, to be proved rigorously below, consider e.g. the se- 
quence AABABBAB.... Here n(l,3) = -1, n(l,7) = -4, n(8,4) = 2, n(3, 3) = 0, and so 
on. Note that if a sequence Q has no AB pair, uab = OVi,j and the matrix G becomes 
identical to the matrix found in [22|. This is expected since the system behaves in this case 
like a single-species TASEP. 

Proof: We have to show first that det G satisfies the master equation. In the absence of 
any nearest neighbor pair, the master equation is 
d ^ 

_pQ|Q t\y; 0) = ^ P^l^ {xi, X2, Xi - 1, Xi+i, Xn; t\y; 0) 
1=1 

_jYpQIQ(x;t|y;0). (29) 

Since all the elements of the matrix G contain the factor e~*, one can write G = e~^G which 
gives 

-f-detG = e-^*^detG-ArdetG. (30) 
dt dt ^ ^ 
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Now the time-derivative in the first term on right-hand side can be written as 



d 
dt 



Gn Gi2 ■■■ Gin 
G21 G22 ••• G2N 



+ 



Gni Gn2 ••• Gnn 
Gil G12 ... Gin 
G21 G22 ■■■ G2N 



Gati Gn2 ■■■ G 



NN 



Gn 


G12 ••• 


Gin 


G21 


G22 ■■■ 


G2N 


Gni 


Gn2 ■■■ 


Gnn 


Gn 


G12 


Gin 


G21 


G22 ■■■ 


G2N 


Gni 


Gn2 ■■■ 


Gnn 


-yi 




{xi - y2 



(31) 



^n(i,N)[Xi — VNj 

where we have used the notation F = e~^F. 

To find the time-derivative of these elements of the i-th row we use 

the identity Fn{x) = —Fn{x) + — 1) from which it follows that 

Fn{x) = Fn{x — 1). Therefore, taking the time derivative generates the row 



Fn(i,i){xi-l-yi),Fri(i,2)ixi-l-y2),:.,Fn(i,N)ixi-l-yN)j- Thus each term on 
the r.h.s. of (l3T]) contributes to one term in the summation in Eq. (129|) . From Eq. (l30!l it 
follows therefore that detG satisfies Eq. fl29|) . 

Now consider the case when one nearest neighbor pair is present: xj = xj-i + 1. If the 
(j — l)-th particles is A and the j-th one is B, then the master equation becomes 



d 



-P^IQ (x; t\y; 0) = J] P'^IQ ^^^^ _ _^ _ ^.^^^ ^|y. g) 



-NpQ\Q{xnX2,...,XN;t\yny;0). 

which gives rise to the boundary condition 

P^^^ (xi, X2, Xj-i, Xj-i, Xj+i, ...,xn] t\-y] 0) = 



(32) 



(33) 



to be satisfied at all times. Since there is an AB pair at (j — 1), j, from the definition of 
n{i,j) it follows that n{i,j — 1) = n{i,j) Vi. Hence in the determinantal representation of 



pQ\Q 



Xi,X2 



, ...,Xj 



-1, xn', tly; 0) the (j — l)-th and j-th row becomes identical 



and the determinant vanishes. If the (j — l)-th particle is B and the j-th one is A then the 
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boundary condition is 

P^^^ {xi, X2, Xj^i, Xj^i, Xj+i, ...,xn; t\y; 0) = (34) 
P^^^ (xi, X2, ...,Xj-i, Xj^i + 1, Xj+i, xn] t\-y] 0) 

which becomes 

Fn{j j -1 - yi) -F„(j_i,2)(a;j„i - 1/2) ••• _ 

Fn(j-l,l)+liXj^l - yi) Fn(j-l^2)+liXj^i - y2) ... 

-Z/i) i^nO-i,2)(a;j-i -^2) -. ^^^^ 

Fn{j-i,i)+i{xj-i + l -yi) -F„(j-i,2)+i(.^i-i + 1 - 2/2) ... 

which can be proved using + = 

If (j — l)-th and j-th particles are AA or BB, then the corresponding boundary condition 
becomes 

P'^l'^ (xi, a;2, Xj-i, Xj+i, xat; t|y; 0) = (36) 
P'^l'^ (xi, X2, Xj_i, Xj_i + 1, Xj+i, Xn; t\y; 0) 

which is easily shown to be satisfied using the identity -F„(x + 1) = Fn{x) — F„_i(x). 

Finally we have to show that det G satisfies the correct initial condition. To this end note 
first the identity Fn{x;t = 0) = 6xfl which follows from the definition of Fn{x;t) in ( |T3l) . 
This means that for xi > yi all the elements of the first column must vanish: Gn = for 
t = which yields detG = 0. For xi = yi one has Gn = 1 and Gn = Vi 7^ 1, since Xj > yi 
for i > 1. Then det G becomes equal to det G^^^ where the matrix G^^'^ is obtained from 
G after omitting the first row and the first column. Next, we assume X2 > I/2 and using the 
above identity, we have G^-J'' = Vi and for X2 = 1/2 G[1^ = 6i^i. It follows therefore that 
detG = 

Sxi,yiSx2,y2 ^^"^^ ^ with G*-^-* being an (A^ — 2) dimensional matrix which results 
from omitting first two rows and columns of G. Repeating the above procedure N times, 
one finds 

N 

detG|i=o = n^-«*- (37) 

i=l 

This completes the proof of the theorem stated above. 
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IV. GEOMETRICAL TREATMENT OF THE BETHE ANSATZ SOLUTION 



In this section, we use the geometrical interpretation of the Bethe Ansatz introduced 



m 



41 



42| to analyze the entangled systems of allowed and forbidden trajectories of different 
species of particles. This approach serves as a heuristic alternative to the explicit rigorous 
solution of the master equation. 

To illustrate this approach let us first consider the single-species TASEP with N = 2. 
The solution of the master equation in this case is 22 1 



Fo(xi - yi)Fo{x2 - y2) - F,{x2 - yi)F_,{x, - y^). (38) 

Instead of the TASEP let us now consider two vicious random walkers (VRW) where two 
particles perform biased random walks on a lattice and if they jump on each other, they 
annihilate. It is easy to verify that the solution of master equation in this case is 

^o(a;i - yi)FQ{x2 - y2) - Fq{x2 - yi)Fo{xi - ys)- (39) 

Notice that in both the TASEP and the VRW, the first term is the product of two Poisson 
processes, which represents two non-interacting random walkers starting from yi,y2 and 
reaching xi,X2, respectively. This term includes all possible trajectories for the pair of 
particles. For the VRW, however, trajectories are not allowed to intersect as they will result 
in annihilation. So one must subtract those forbidden trajectories where there exists at least 
one intersection. This is obtained by considering two non-interacting random walkers again, 
one starting from yi, reaching X2 and the other starting from y2 and reaching Xi. Since 
Xi < X2 and yi < y2 their trajectories must cross at least once and therefore (l39l) yields the 
correct transition probability. This observation follows also directly from the free-fermion 



representation of this process, see [2, l43l446l| for a detailed treatment. 

For a TASEP, however, if the particles attempt to jump onto each other, they are not 
annihilated but are reflected off each other (see Fig[T]), i.e., the jump attempt is unsuccessful. 
This means out of all the trajectories forbidden for the VRW we allow those paths for the 
TASEP where particles collide and get reflected off. As follows from the comparison of Eqs. 
( l38l) and ( l39i) . the weight of this type of trajectories is Fo{x2 — yi)Fo{xi — 1/2) — Fi{x2 — 
?/i)F_i(xi - y2). 

Now consider one A particle and one B particle. First consider the probability that the 
initial sequence {AB} remains unchanged at time t. If the two particles jump on each other, 
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A 



A 



y 




A 



A 



FIG. 1: Mean trajectories in the single-species TASEP. When trajectories of two particles collide 
they reflect off each other. Here the horizontal line at the top shows the configuration of the system 
at time t and that at the bottom shows the initial configuration. 

they exchange positions. Once the exchange has taken place, the particles are not allowed 
to go back and the sequence remains BA thereafter. This means that if we want to compute 
the probability that even at time t the sequence continues to be {AB}, then we must not 
allow the particles to jump on each other. All those trajectories where a jump attempt has 
been made, must be forbidden. In other words, the two particles must behave like a pair of 
VRW's. This explains the form of the solution in ( fT5l) . 

For the case when the A and B particles actually interchange and the sequence becomes 
BA, one must consider only those paths where there exists a crossing, as shown in Fig |2j 
Note that this diagram can be obtained from Fig [T] after 'recoloring' the lines. Even for 
the single-species TASEP, where particles are not allowed to interchange positions, one can 
reinterpret the diagram[T]as a process where two particles actually interchange positions and 
then renamed after the interchange (corresponds to the exchange of the colors of the lines). 
Hence Figures [1] and [2] are essentially identical and the weight of the allowed trajectories in 
this case is 



as derived above. Again, our rigorous solution in ( lT5l) conforms to this argument. 

Thus we have developed an equivalent diagrammatic approach to interpret the Bethe 
Ansatz solution presented in section [Tll This approach can easily be generalized for higher 
values of A^. There are three different possibilities when two particles interact — when an 
A interacts with a i? on its right, they can behave like VRW or they can exchange. For 
any other collision, the particles are reflected off each other (in Fig [T] the trajectory starting 



Fo{x2 - yi)Fo{xi - 1/2) - Fi{x2 - ?/i)F_i(xi - 1/2) 



(40) 
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A B 

FIG. 2: Diagram showing interchange of an A and B particle 

from yi is reflected backward while the other starting from 1/2 is reflected forward). For a 
set of particles, there are A^! number of possible diagrams obtained from connecting all 
{yi, y2, 2/Ar} with all {a;i, X2, xn}- Clearly, each diagram corresponds to one particular 
permutation for the Bethe momentum variables {pi,P2, ■■■,Pn} with the conjugate position 
variables {xi,X2, ...,xn}- In a diagram, if yj is joined to Xi, then the corresponding term 
in the solution must have an F function with argument (x, — yj). To obtain the index 
(subscript) of this F function, count the number of intersections on the line joining Xi and 
yj. Each forward (backward) reflection corresponds to an index -1(+1), each VRW collision 
corresponds to an index and each AB interchange is represented as difference of two terms, 
as in Eq. (HOl) . For illustration, we provide an example below. 

In Fig [3] the intersection marked V is between an AB pair and since they do not inter- 
change they must behave as VRW. So collision V yields an index in the F functions with 
arguments (xi — y^) and {x^ — 7/2)- The collision R is between two A particles and hence of 
ordinary reflecting type, where the line {x2 — yi) suffers a backward reflection and {xi — ys) 
has a forward reflection. Final weight of the diagram is Fq^x^ — y2)F^i{xi — y3)Fi{x2 — yi). 

Consider another example. In Fig H] in the collision I the A and B interchange and in 
collision R the line {x2 — yi) has a backward reflection and {xi — y^) has forward reflection. 
These taken together yields Fi{x2-yi) [F^iixi - y3)Fo{x3 - 1/2) - ^^-1(3:3 - y2)Fo{xi - y^)]. 
Note that the part in the square bracket has the same form as in Eq. fHOj) but because 
[x\ — 7/3) suffers a forward reflection with (x2 — yi), its index is lowered by unity. 

Note that not all the A^! diagrams are always allowed. In the above example, when an 
initial sequence AAB changed to ABA, then not all diagrams are possible: Since we want 
the B particle starting at 2/3 to interact with the other A particles, we cannot have any 
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FIG. 3: The sequence AAB remains same in time and the AB pair must behave hke VRW. 




A A B 

FIG. 4: Sequence AAB changes to ABA and in the colhsion marked / the AB pair actually 
interchange and collision R is of reflecting type 

diagram where is joined to ?/3, meaning no intersection lies on the (0:3 — y^) line. In Bethe 
Ansatz language this means that all permutations where P3X3 appears are absent from the 
solution. 

Thus these diagrams serve as an alternative way to interpret the terms in the Bethe 
Ansatz solution. This approach also offers us an insight into the theorem discussed in 
section lllli We can now understand why the determinant in [22| for the single species 
case gets modified into det G for the two-species case obtained above. When only one type 
of particles is present, the intersections in their trajectories result in only reflecting type 
collisions and the corresponding index of the F function is either +1 or —1. It can be easily 
seen that the line (xj — Uj) picks up an index {i — j) for the F function, as a result of all 
these reflections. But when both A and B particles are present, and if we are interested 
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in the case when their sequence do not change, then some of these coUisions have to be of 
'vicious' type which corresponds to an index 0. For the hne (xj — Uj) the number of such 
'vicious' coUisions is uab-, the number of AB pairs between the z-th and j-th particle. As a 
result, the corresponding index of the F function is not {i — j) but sgn{i — — j \ — uab)- 
This explains the structure of the matrix G in section lllli 

V. FURTHER DISCUSSION 

In this paper we have used the Bethe Ansatz to solve the master equation of the TASEP 
with first and second class particles for arbitrary initial conditions. We have also obtained a 
compact determinantal representation of the solution for the case when the initial sequence 
of particles remains unchanged. An equivalent geometrical approach developed in section 
IIVI provides us with insight into the form of the determinant. We have seen that when the 
sequence does not change with time, then all the AB pairs in the sequence must behave like 
VRWs. 

Interestingly the last conclusion remains valid even when more than two species of parti- 
cles are present in the system. For any initial sequence Q, consisting of any arbitrary number 
of species (classes) of particles, if we are interested in the probability that the sequence does 
not change until time t, then all the pairs, where a higher class particle has a lower class 
particle to its right, must behave like VRWs. All the trajectories where they jump onto each 
other, must therefore be forbidden because once they interchange positions, they cannot go 
back again. Generalizing the geometrical approach of section HVt it is easy to see that even in 
this case the diagrams contain two different types of intersection points: one of vicious type 
(corresponding to an index 0) and the other of reflecting type (corresponding to an index 
±1). For any given line (xj — Vj), there will be an F function as before, with the argument 
{xi — yj) and index n(i,j) = sgn{i — — j\ —n^) where Uy denotes the number of 'vicious 
pairs' (where a higher class particle has a lower class particle on its right) between i-th and 
j-th member of the sequence. Thus using this approach one can write down a determinantal 
representation for the solution of the master equation for a general multi-class TASEP where 
the sequence of various species of particles remains unchanged. 

In the special case of a sequence where the first particle from the left is of class 1, the 
second particle is of class 2, the third particle is of class 3, ...,the iV-th particle is of class 
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A^, it follows from the above argument that all the pairs in this sequence behave like VRW's 
if we consider the probability that the sequence is preserved until time t. This probability 
is nothing but det where is a matrix with (zj)-th element F^ij = Fo{xi — yj). If 
on the other hand the initial sequence is reversed when the first class particle becomes the 
rightmost particle and the A^-th class particle becomes the leftmost particle, then there are 
no vicious pairs left in the system and the solution for the single species TASEP remains 
valid. Thus depending on the labels of the particles in the initial sequence, one can have 
TASEP solution or VRW solution. 

In this context it is interesting to discuss another recent result on the crossover from VRW 
to TASEP. In a model for semi-vicious walkers was introduced that interpolates between 
VRW and TASEP, having the two as limiting cases. In this model two particles annihilate 
if they jump on each other, but this jump attempt to an occupied site takes place with a 
reduced probability. Note that for VRW's the probability to jump on an occupied site is same 
as that for an empty site, while for the TASEP jumping on an occupied site is prohibited. 
Thus by varying the probability to jump on an occupied site one can go from one limiting 
case to the other and the model for semi-vicious walkers serves an intermediate between 
these two limits. It was shown that the survival probability for A^ particles can be described 



by a scaling form that characterizes the transition from VRW to TASEP [47| . In contrast we 
consider here a multi-species exclusion process and have argued in the previous paragraph 
that if each of the A^ particles belongs to a distinct class, then it is possible to have a VRW 
or a TASEP solution depending on the label of the particles in the initial sequence. For 
a general sequence, under the condition that until time t the sequence remains unchanged, 
some of the particles behave like TASEP and some behave like VRW. Thus both these 
aspects are simultaneously present in our model. In the case when the number of vicious 
interfaces Uy = N — 1, all the particles behave like VRW and the solution of the master 
equation becomes det F^, while for = all the particles behave as TASEP particles. For 
intermediate values of both these limiting behaviors are present and some pairs behave 
like TASEP and some like VRW and the determinantal solution of the master equation is 
modified accordingly. 

For the general two-species problem, i.e. when the initial sequence of species changes 
because of interchange of A and B particles, it is an interesting question whether the Bethe 
ansatz solution presented in this work allows for a compact representation in terms of de- 
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terminants. At this point we do not have a complete answer to this question. We have been 
able to obtain compact representations for certain special cases. For example, if the final 
sequence Q' results from a single AB interchange from the initial sequence Q, such that 
the number of AB pairs in Q' is less than that in Q, then we have been able to prove that 
the probability P^'l<3 = P«'IQ' - pQ\Q can be written difference of two determinants. 
The existence of similar compact representation for any number of interchanges and for the 
general multi-class TASEP is an open problem for future investigation. 
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Appendix A: Verification of Yang-Baxter Criterion 

In this appendix, we discuss the validity of Yang-Baxter equation for = 3 case. Con- 
sider the r.h.s. of Eq. ( IT71) . The matrices Sf2(Pj,PA:) and Ti^^^ipj^pk) represents elementary 
permutations of the momenta variables. For example, T,l2{Pj,Pk) represents the interchange 
of the momentum variables pj and pk which are coupled to Xi and X2, respectively. Each 
term in Eq. ( IT7|) corresponds to one particular permutation V of the momenta and the 
product of the S matrices represents decomposition of V into a series of elementary per- 
mutations. This decomposition need not be unique for all V. Consider, for example, the 
last term in Eq. ( TTTl) . where the momenta variables coupled to xi,X2,X3, are, respectively, 
P3,P2,Pi- Starting from a permutation Pi,P2,P3 (as in the first term on the r.h.s. of Eq. 
( !T7|l . to reach the arrangement P3,P2,Pi, one can first interchange (^25^3)5 then (^1,^3) and 
again (^1,^2)- Accordingly, the S-product is T^^^^pi, p2)T,l2{pi, P3)^23iP'2^ Ps) ■ But the same 
permutation can also be reached by first interchanging (^1,^2), then (^1,^3) and finally 
(^2,^3)- In that case, the corresponding matrix will be 2 (p2,P3)S23(pi,P3)Sf 2(^1,^2)- In 
order for the system to be integrable, one must satisfy the Yang-Baxter criterion 



^12 (P2 ,^3)^23 (Pl , P3) S12 (Pl , P2 ) = (pi , P2 ) 2^2 (Pl ,^3)^23(^2,^3)- 



(Al) 
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(We remark that the shghtly different Yang-Baxter equation for the scattering matrix S 
is obtained from this croterion through the relation S = HP where P is the permutation 
operator). Now, using the definitions T,l2{Pj,Pk) = ^{Pj,Pk) ® I2 and T,2^{pj,pk) = I2 ® 
T.{pj,pk), we have 
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where Sjk = -(1 - e'P^)/{l - e'^^) and Tjk = -{e'P" - e'P^)/{l - e'^^). Using these forms it 
is straightforward to verify Eq. flAip . 

We remark that the nested Bethe ansatz can be also be used for treating second-class 
particles in the partially asymmetric simple exclusion process (PASEP) where particles jump 



to the right or left with non-zero rates 



However, even for the single-species PASEP 



there is no known determinantal representation of the transition probabilities. 
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